Simulation of the Undiseased Human Cardiac Ventricular Action Potential: Model Formulation and Experimental Validation

Cellular electrophysiology experiments, important for understanding cardiac arrhythmia mechanisms, are usually performed with channels expressed in non myocytes, or with non-human myocytes. Differences between cell types and species affect results. Thus, an accurate model for the undiseased human ventricular action potential (AP) which reproduces a broad range of physiological behaviors is needed. Such a model requires extensive experimental data, but essential elements have been unavailable. Here, we develop a human ventricular AP model using new undiseased human ventricular data: Ca2+ versus voltage dependent inactivation of L-type Ca2+ current (ICaL); kinetics for the transient outward, rapid delayed rectifier (IKr), Na+/Ca2+ exchange (INaCa), and inward rectifier currents; AP recordings at all physiological cycle lengths; and rate dependence and restitution of AP duration (APD) with and without a variety of specific channel blockers. Simulated APs reproduced the experimental AP morphology, APD rate dependence, and restitution. Using undiseased human mRNA and protein data, models for different transmural cell types were developed. Experiments for rate dependence of Ca2+ (including peak and decay) and intracellular sodium ([Na+]i) in undiseased human myocytes were quantitatively reproduced by the model. Early afterdepolarizations were induced by IKr block during slow pacing, and AP and Ca2+ alternans appeared at rates >200 bpm, as observed in the nonfailing human ventricle. Ca2+/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca2+ cycling. INaCa linked Ca2+ alternation to AP alternans. CaMK suppression or SERCA upregulation eliminated alternans. Steady state APD rate dependence was caused primarily by changes in [Na+]i, via its modulation of the electrogenic Na+/K+ ATPase current. At fast pacing rates, late Na+ current and ICaL were also contributors. APD shortening during restitution was primarily dependent on reduced late Na+ and ICaL currents due to inactivation at short diastolic intervals, with additional contribution from elevated IKr due to incomplete deactivation

Mechanisms of persistent atrial fibrillation and recurrences within 12 months post-ablation: Non-invasive mapping with electrocardiographic imaging

INTRODUCTION: Catheter ablation of persistent AF has not been consistently successful in terminating AF or preventing arrhythmia recurrences. Non-invasive Electrocardiographic Imaging (ECGI) can help to understand recurrences by mapping the mechanisms of pre-ablation AF and comparing them with the patterns of recurrent arrhythmias in the same patient. METHODS: Seventeen persistent AF patients underwent ECGI before their first catheter ablation. Time-domain activation maps and phase progression maps were obtained on the bi-atrial epicardium. Location of arrhythmogenic drivers were annotated on the bi-atrial anatomy. Activation and phase movies were examined to understand the wavefront dynamics during AF. Eight patients recurred within 12 months of ablation and underwent a follow-up ECGI. Driver locations and movies were compared for pre- and post-ablation AF. RESULTS: A total of 243 focal drivers were mapped during pre-ablation AF. 62% of the drivers were mapped in the left atrium (LA). The pulmonary vein region harbored most of the drivers (43%). 35% of the drivers were mapped in the right atrium (RA). 59% (10/17) and 53% (9/17) of patients had repetitive sources in the left pulmonary veins (LPV) and left atrial appendage (LAA), and the lower half of RA, respectively. All patients had focal drivers. 29% (5/17) of patients had macro-reentry waves. 24% (4/17) of patients had rotors. Activation patterns during persistent AF varied from single macro-reentry to complex activity with multiple simultaneous wavefronts in both atria, resulting in frequent wave collisions. A total of 76 focal driver activities were mapped in 7/8 patients during recurrence. 59% of the post-ablation AF drivers were mapped in the LA. The pulmonary vein region harbored 50% of total drivers. 39% of sources were mapped in the RA. AF complexity remained similar post-ablation. 58% (44/76) of pre-ablation sources persisted during recurrence. 38% (3/8) of patients had macro-reentry and one patient had rotors. CONCLUSION: ECGI provides patient-specific information on mechanisms of persistent AF and recurrent arrhythmia. More than half pre-ablation sources repeated during post-ablation recurrence. This study provides direct evidence for drivers that persist days and months after the ablation procedure. Patient-tailored bi-atrial ablation is needed to successfully target persistent AF and prevent recurrence. ECGI can potentially predict recurrence and assist in choice of therapy

Noninvasive electrocardiogram imaging of substrate and intramural ventricular tachycardia in infarcted hearts

AbstractObjectivesThe goal of this study was to experimentally evaluate a novel noninvasive electrocardiographic imaging modality during intramural reentrant ventricular tachycardia (VT).BackgroundMyocardial infarction and subsequent remodeling produce abnormal electrophysiologic substrates capable of initiating and maintaining reentrant arrhythmias. Existing noninvasive electrocardiographic methods cannot characterize abnormal electrophysiologic substrates in the heart or the details of associated arrhythmias. A noninvasive method with such capabilities is needed to identify patients at risk of arrhythmias and to guide and evaluate therapy.MethodsA dog heart with a four-day-old infarction was suspended in a human shaped torso-tank. Measured body surface potentials were used to noninvasively compute epicardial potentials, electrograms and isochrones. Accuracy of reconstruction was evaluated by direct comparison to measured data. Reconstructions were performed during right atrial pacing and nine cycles of VT.ResultsNoninvasively reconstructed potential maps, electrograms and isochrones identified: 1) the location of electrophysiologically abnormal infarct substrate; 2) the epicardial activation sequences during the VTs; 3) the locations of epicardial breakthrough sites; and 4) electrophysiologic evidence for activation of the Purkinje system and septum during the reentrant beats.ConclusionsElectrocardiographic imaging can noninvasively reconstruct electrophysiologic information on the epicardium during VT with intramural reentry, provide information about the location of the intramural components of reentry and image abnormal electrophysiologic substrates associated with infarction

Characterization of double potentials in a functionally determined reentrant circuit Multiplexing studies during interruption of atrial flutter in the canine pericarditis model

AbstractObjectives. We tested the hypothesis that double potentials recorded during atrial flutter in a functionally determined reentrant circuit reflect activation of the reentrant wave front around an area of functional conduction block.Background. The center of the atrial flutter reentrant circuit in the sterile pericarditis canine model is characterized by double potentials.Methods. We studied 11 episodes of atrial flutter in eight dogs during interruption of atrial flutter while pacing the atria. A multielectrode mapping system was used to record simultaneously from 190 electrodes on the right atrium (location of reentry).Results. Interruption of atrial flutter occurred when the orthodromic wave front from the pacing impulse blocked in an area of slow conduction in the reentrant circuit. The response of the double potential with interruption of atrial flutter depended on the location of the recording site relative to this area of block. Two types of response were seen. When the double potential was recorded orthodromically distal to this area of block, interruption of atrial flutter was associated with disappearance of the second deflection, and continued pacing after interruption of atrial flutter was not associated with reappearance of the second potential. When the double potential was recorded at a site orthodromically proximal to the area of block, interruption of atrial flutter was not associated with disappearance of the second potential, and when rapid atrial pacing was continued, the double potential remained despite disappearance of the atrial flutter reentrant circuit.Conclusions. Double potentials represent functional conduction block in the center of the reentrant circuit, with each deflection of the double potential reflecting activation on either side of the area of functional block. The data also demonstrate that double potentials are not limited to a reentrant circuit, as they were recorded on either side of an area of block in the absence of such a circuit

Determinants of Beat-to-Beat Variability of Repolarization Duration in the Canine Ventricular Myocyte: A Computational Analysis

Beat-to-beat variability of repolarization duration (BVR) is an intrinsic characteristic of cardiac function and a better marker of proarrhythmia than repolarization prolongation alone. The ionic mechanisms underlying baseline BVR in physiological conditions, its rate dependence, and the factors contributing to increased BVR in pathologies remain incompletely understood. Here, we employed computer modeling to provide novel insights into the subcellular mechanisms of BVR under physiological conditions and during simulated drug-induced repolarization prolongation, mimicking long-QT syndromes type 1, 2, and 3. We developed stochastic implementations of 13 major ionic currents and fluxes in a model of canine ventricular-myocyte electrophysiology. Combined stochastic gating of these components resulted in short- and long-term variability, consistent with experimental data from isolated canine ventricular myocytes. The model indicated that the magnitude of stochastic fluctuations is rate dependent due to the rate dependence of action-potential (AP) duration (APD). This process (the “active” component) and the intrinsic nonlinear relationship between membrane current and APD (“intrinsic component”) contribute to the rate dependence of BVR. We identified a major role in physiological BVR for stochastic gating of the persistent Na+ current (INa) and rapidly activating delayed-rectifier K+ current (IKr). Inhibition of IKr or augmentation of INa significantly increased BVR, whereas subsequent β-adrenergic receptor stimulation reduced it, similar to experimental findings in isolated myocytes. In contrast, β-adrenergic stimulation increased BVR in simulated long-QT syndrome type 1. In addition to stochastic channel gating, AP morphology, APD, and beat-to-beat variations in Ca2+ were found to modulate single-cell BVR. Cell-to-cell coupling decreased BVR and this was more pronounced when a model cell with increased BVR was coupled to a model cell with normal BVR. In conclusion, our results provide new insights into the ionic mechanisms underlying BVR and suggest that BVR reflects multiple potentially proarrhythmic parameters, including increased ion-channel stochasticity, prolonged APD, and abnormal Ca2+ handling

State-dependent electrostatic interactions of S4 arginines with E1 in S2 during Kv7.1 activation

The voltage-sensing domain of voltage-gated channels is comprised of four transmembrane helices (S1–S4), with conserved positively charged residues in S4 moving across the membrane in response to changes in transmembrane voltage. Although it has been shown that positive charges in S4 interact with negative countercharges in S2 and S3 to facilitate protein maturation, how these electrostatic interactions participate in channel gating remains unclear. We studied a mutation in Kv7.1 (also known as KCNQ1 or KvLQT1) channels associated with long QT syndrome (E1K in S2) and found that reversal of the charge at E1 eliminates macroscopic current without inhibiting protein trafficking to the membrane. Pairing E1R with individual charge reversal mutations of arginines in S4 (R1–R4) can restore current, demonstrating that R1–R4 interact with E1. After mutating E1 to cysteine, we probed E1C with charged methanethiosulfonate (MTS) reagents. MTS reagents could not modify E1C in the absence of KCNE1. With KCNE1, (2-sulfonatoethyl) MTS (MTSES)− could modify E1C, but [2-(trimethylammonium)ethyl] MTS (MTSET)+ could not, confirming the presence of a positively charged environment around E1C that allows approach by MTSES− but repels MTSET+. We could change the local electrostatic environment of E1C by making charge reversal and/or neutralization mutations of R1 and R4, such that MTSET+ modified these constructs depending on activation states of the voltage sensor. Our results confirm the interaction between E1 and the fourth arginine in S4 (R4) predicted from open-state crystal structures of Kv channels and reveal an E1–R1 interaction in the resting state. Thus, E1 engages in electrostatic interactions with arginines in S4 sequentially during the gating movement of S4. These electrostatic interactions contribute energetically to voltage-dependent gating and are important in setting the limits for S4 movement

Study protocol: MyoFit46-the cardiac sub-study of the MRC National Survey of Health and Development

BACKGROUND: The life course accumulation of overt and subclinical myocardial dysfunction contributes to older age mortality, frailty, disability and loss of independence. The Medical Research Council National Survey of Health and Development (NSHD) is the world's longest running continued surveillance birth cohort providing a unique opportunity to understand life course determinants of myocardial dysfunction as part of MyoFit46-the cardiac sub-study of the NSHD. METHODS: We aim to recruit 550 NSHD participants of approximately 75 years+ to undertake high-density surface electrocardiographic imaging (ECGI) and stress perfusion cardiovascular magnetic resonance (CMR). Through comprehensive myocardial tissue characterization and 4-dimensional flow we hope to better understand the burden of clinical and subclinical cardiovascular disease. Supercomputers will be used to combine the multi-scale ECGI and CMR datasets per participant. Rarely available, prospectively collected whole-of-life data on exposures, traditional risk factors and multimorbidity will be studied to identify risk trajectories, critical change periods, mediators and cumulative impacts on the myocardium. DISCUSSION: By combining well curated, prospectively acquired longitudinal data of the NSHD with novel CMR-ECGI data and sharing these results and associated pipelines with the CMR community, MyoFit46 seeks to transform our understanding of how early, mid and later-life risk factor trajectories interact to determine the state of cardiovascular health in older age. TRIAL REGISTRATION: Prospectively registered on ClinicalTrials.gov with trial ID: 19/LO/1774 Multimorbidity Life-Course Approach to Myocardial Health- A Cardiac Sub-Study of the MCRC National Survey of Health and Development (NSHD)